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Abstract 

This  report  deals  with  the  manifestation  of  downslope  converted  continental  shelf  noise 
sources  in  the  matched-field  processor  (MFP)  ambiguity  surface.  Of  interest  here  is  how 
noise  sources  outside  the  range  of  interest  leak  into  the  ambiguity  surface  through  the 
sidelobe  structure  of  the  effective  MFP  beam  (or  cell)  pattern.  The  focus  is  on  the  appear¬ 
ance  of  the  matched-field  processing  ambiguity  surface  due  to  continental  shelf  noise,  pri¬ 
marily  ships,  that  gets  coupled  in  the  deep  sound  channel.  This  study  is  carried  out  with 
simulations  using  the  parabolic  equation  (PE)  model.  PE  is  used  to  simulate  the  pressure 
field  across  a  receiving  vertical  line  array  in  deep  water  (4000  m),  1750  km  from  the  con¬ 
tinental  shelf.  A  shallow,  low  frequency  (18  Hz)  source  is  assumed  deployed  at  various 
locations  on  the  continental  slope.  The  MFP  ambiguity  surfaces  are  calculated  on  a  search 
window  extending  in  range  across  500  km  and  in  depth  down  to  1000  m,  and  obtained 
with  a  1000  m  or  a  4000  m  long  simulated  receiving  array.  Shallow  sources  on  the  slope 
show  up  in  the  MFP  ambiguity  surface  as  deep  peaks  when  downslope  conversion  takes 
place.  Shallow  sources  at  locations  on  the  slope  or  in  deep  water  adjacent  to  the  slope 
where  downslope  conversion  does  not  take  place,  appear  in  the  MFP  ambiguity  surface  as 


shallow  sources. 
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1.  Introduction 

This  report  deals  with  the  manifestation  of  downslope  converted  continental  shelf 
noise  sources  in  the  matched-field  processor  (MFP)  ambiguity  surface. 

Downslope  conversion  of  acoustic  energy  has  been  suggested  as  one  mechanism 
by  which  noise  radiated  by  ships  or  storms  gets  coupled  into  the  sound  channel.  Various 
experiments  have  been  designed  to  better  understand  downslope  conversion,  such  as  the 
one  conducted  by  the  Marine  Physical  Laboratory  (MPL)  in  July  1989  in  the  region  of  the 
continental  margin  NW  of  Pt  Conception,  California  [1]. 

Measurements  of  propagation  loss  show  that  bathymetry  does  strongly  influence 
the  propagation  of  sound  from  a  shallow  source  and  downslope  enhancements  of  up  to  10 
dB  were  observed  at  a  deep  receiver  for  sources  near  the  edge  of  the  continental  shelf  [2- 
4].  Downslope  conversion  also  provides  an  explanation  for  the  vertical  directionality  of 
low  frequency  ambient  noise  in  the  deep  ocean  [5-6].  The  grazing  angle  of  acoustic  rays 
propagating  upon  a  sloping  bottom  are  reduced  by  two  times  the  slope  angle  upon  refec¬ 
tion.  Hence,  a  shallow  source  on  the  continental  shelf  produces  arrivals  in  deep  water  with 
low  angles  [7].  At  long  distance,  this  downslope  converted  energy  is  seen  as  a  large  contri¬ 
bution  to  the  ambient  noise  level  arriving  near  the  horizontal  [5,6,8]. 

To  complement  observations  of  downslope  conversion  with  transmission  loss  and 
conventional  beamforming  results,  this  report  deals  with  the  manifestation  of  downslope 
converted  continental  shelf  noise  sources  in  matched-field  processing  (MFP)  ambiguity  sur¬ 
faces.  Matched-field  processing  is  a  generalized  beamforming  technique  which  correlates 
the  field  observed  by  an  array  of  receivers  with  the  predicted  field  due  to  a  source  in  a 
particular  range  and  depth  cell.  Of  interest  here  is  the  ambiguity  surface  contribution  due 
to  a  noise  source  outside  the  range-depth  search  window,  in  the  continental  shelf  region, 
which  leaks  through  the  sidelobes  of  the  MFP  "beam"  steered  on  a  particular  hypothesized 
source  location  in  deep  water. 
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The  work  performed  and  presented  in  this  document  deals  with  point  source  simu¬ 
lations  only.  The  range  dependent  Parabolic  Equation  (PE)  model  [9-11]  is  used  to  simulate 
pressure  fields  observed  at  a  vertical  line  array  due  to  low  frequency  sources  placed  at 
various  positions  on  the  continental  slope  as  in  [8].  PE  is  also  used  to  compute  the  set  of 
MFP  replica  vectors,  i.e.  the  predicted  pressure  across  the  array  for  point  sources  in  a 
range-depth  search  window.  The  geometry  and  environmental  information  used  in  these 
simulations  are  similar  to  those  encountered  during  the  MPL  Downslope  Conversion  Exper¬ 
iment  [1]. 

2.  Description  of  the  Simulations 

The  overall  propagation  geometry  is  given  in  Fig.  1.  The  array  of  receivers  is 
assumed  deployed  at  a  fixed  location  in  deep  water  (nominal  depth  of  4000  m).  The  sub¬ 
stantial  amount  of  environmental  data,  collected  during  the  actual  experiment  [12]  (such  as 
CTDs,  XBTs,  range  varying  temperature  field  obtained  with  AXBTs),  allows  the  use  of  a  real¬ 
istic  range  dependent  sound  speed  structure  from  the  array  to  the  continental  slope  in  this 
simulation  study.  Together  with  the  bathymetry,  the  sound  speed  field  is  plotted  in  Fig.  2 
for  the  deep  water  and  in  Fig.  3  for  the  slope  area. 

The  first  deep  water  environmental  sector  starts  from  the  receivers  and  is  30  km 
wide.  It  is  followed  by  27  other  deep  water  sectors  (which  are  60  km  wide),  to  the  con¬ 
tinental  slope.  The  deep  water  sectors  have  a  nominal  water  depth  of  4000  m.  The  transi¬ 
tion  from  deep  water  to  the  continental  shelf  begins  with  a  gentle  0.4  deg  slope  for  57.5 
km,  followed  by  a  sharp  rise  over  8.8  km  from  a  depth  of  3600  m  to  a  depth  of  1800  m. 
This  steep  slope  is  followed  by  a  small  2.8  deg  slope  and  a  shelf  with  a  depth  of  480  m. 

The  sound  speed  field  is  strongly  range  dependent  and  Fig.  2  clearly  shows  the 
shoaling  of  the  sound  axis  towards  the  continental  shelf.  The  geoacoustic  model  for  the 
28  deep  water  sectors  corresponds  to  a  fairly  reflecting  bottom  where  the  350  m  thick  sed¬ 
iment  layer  has  a  sound  speed  increasing  from  1523.6  m/s,  the  bottom  water  sound  speed. 
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to  3500  m/s.  The  density  is  2  g/cm^  and  the  attenuation  is  0.113  dB/wavelength.  These 
parameters  attempt  to  simulate  an  abyssal  plain  with  a  very  thin  sediment  layer  on  top  of 
hard  basalt.  The  two  wavelength  thick  PE  subbottom  has  the  same  properties  as  the  sedi¬ 
ments  with  a  sound  speed  equal  to  3500  m/s. 

The  bottom  sound  speed  profiles  in  the  continental  slope  area  are  given  in  Fig.  4. 
The  shallow  water  sector  sediment  layer  has  a  density  of  1.5  g/cm^,  and  an  attenuation  of 
0.113  dB/wavelength.  The  steep  slope  sector  starting  at  -33.7  km  in  Fig.  3  is  assumed  to 
have  no  sediment  deposits  and  the  sound  speed  is  constant  equal  to  2500  m/s,  the  density 
is  2  g/cm^,  the  attenuation  0.113  dB/wavelength.  The  gentle  slope  sector  starting  at  -42.5 
km  in  Fig.  3  has  a  density  equal  to  1.5  g/cm^  and  an  attenuation  of  0.1 13  dB/wavelength. 

An  18  Hz  source  is  assumed  deployed  at  shallow  depth,  at  four  different  locations 
around  the  continental  slope.  The  source  is  first  on  the  continental  shelf  in  shallow  water 
(480  m),  1750  km  from  the  receivers  (at  range  0  in  Fig.  3).  It  is  then  moved  in  the  middle  of 
the  small  2.8  deg.  slope,  1729.6  km  from  the  array  (at  range  -20.4  km  in  Fig.  3).  The  third 
station  corresponds  to  the  middle  of  the  steep  11.5  deg.  slope  sector,  1711.9  km  from  the 
array  (in  Fig.  3  at  range  -38.1  km).  Finally,  the  source  is  set  at  the  beginning  of  the  gentle 
0.4  deg  slope  sector,  1707.5  km  from  the  fixed  receivers  (in  Fig.  3  at  range  -42.5  km). 

The  modified  wide  angle  PE  model  [10]  is  run  with  a  1024  point  transform  (the 
sampling  in  depth  is  on  the  order  of  6.17  m).  The  range  step  is  100  m,  close  to  a 
wavelength  at  18  Hz.  The  starting  field  of  a  source  set  at  a  depth  of  45  m  (approximately 
half  a  wavelength  at  18  Hz)  is  used  to  simulate  a  ship  and  is  calculated  with  the  normal 
mode  model  KRAKEN  [13].  The  receiving  array  is  assumed  fixed  in  range  with  respect  to 
the  environment,  while  the  source  is  moved  from  the  shallow  water  station  to  the  deep 
water.  Two  vertical  line  arrays  are  considered  : 

(1)  an  aperture  on  the  order  of  1000  m,  deployed  from  400  to  1390  m  with  80  hydro¬ 
phones  equally  spaced  every  12.353  m. 
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(2)  a  4000  m  aperture  covering  the  full  water  column  with  323  sensors  spaced  every 

12.353  m. 

3.  Modeling  Results 

In  this  section,  the  simulated  data  are  studied  and  their  physical  validity  verified. 
This  is  carried  out  by  looking  at  transmission  loss  and  field  plots,  as  well  as  conventional 
beamforming  results. 

First,  the  existence  of  the  downslope  enhancement  effect  is  verified.  Fig.  5  gives 
the  transmission  loss  curves  for  a  1000  m  receiver,  as  a  function  of  range  for  the  four  sta¬ 
tions.  The  range  origin  in  Fig.  5  corresponds  to  the  source  in  shallow  water.  The  receiving 
array  is  1750  km  away.  The  three  other  simulated  stations  have  increasing  ranges  equal  to 
20.4  km  (small  2.8  deg  slope),  38.1  km  (the  steep  11.5  deg  slope)  and  42.5  km  (the  begin¬ 
ning  of  the  gentle  slope)  from  the  shallow  water  station.  The  smallest  transmission  loss  at 
the  array  range  (top  curve  in  Fig.  5),  which  is  on  the  order  of  100  dB,  corresponds  to  the 
source  in  the  middle  of  the  small  (2.8  deg)  slope  sector.  The  next  smallest  transmission 
loss  at  the  array  range  (1750  km),  on  the  order  of  105  dB,  corresponds  to  the  source  in 
shallow  water.  The  largest  transmission  losses  (107-108  dB)  occur  when  the  source  is 
above  the  steep  slope  or  at  the  beginning  of  the  gentle  slope.  Downslope  enhancement  is 
therefore  observed  with  levels  on  the  order  of  8  dB,  consistent  with  enhancement  reported 
in  the  literature  of  up  10  dB  [2-4]. 

These  preliminary  results  indicate  that  slope  interaction  happens  when  the  source 
is  either  in  shallow  water  or  above  the  small  2.8  deg.  slope.  A  convenient  way  to  verify 
this  conclusion  is  to  look  at  field  plots  for  the  four  different  simulated  transmission 
experiments.  Field  plots  extending  in  depth  to  4350  m  (hence  covering  the  whole  water 
column  and  the  350  m  thick  sediment  layer  as  well)  and  extending  in  range  up  to  150  km 
from  the  source  are  given  in  Figs.  6,  8,  10  and  12  for  the  four  stations  while  the 
corresponding  bathymetry  is  given  in  Figs.  7,  9,  11,  13.  Slope  interaction  and  energy 
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conversion  from  steep  to  near  horizontal  paths  are  clearly  observed  for  the  source  in  shal¬ 
low  water  in  Fig.  6  and  above  the  small  slope  area  in  Fig.  8.  Higher  levels  are  observed  in 
the  latter  case.  On  the  other  hand,  there  is  no  slope  interaction  when  the  source  is  above 
the  steep  slope  or  at  the  beginning  of  the  gentle  slope.  The  corresponding  field  plots  in 
Figs.  10  and  12  show  steep  bottom  interacting-surface  reflected  paths  occurring  every  50 
km,  forming  the  usual  convergence  zone  pattern  observed  for  shallow  sources  in  deep 
water. 

The  vertical  directionality  of  the  field  received  at  the  array  also  can  be  used  to 
observe  the  effect  of  downslope  conversion.  The  results  of  conventional  FFT  beamforming 
on  the  1000  m  aperture,  given  in  Fig.  14,  clearly  shows  downslope  conversion  when  the 
source  is  either  in  shallow  water  or  above  the  small  2.8  deg  slope.  Significant  power  levels 
are  observed  at  low  angles  (within  ±10  deg.).  There  is  no  power  at  low  arrival  angle 
observed  in  the  vertical  directionality  when  the  source  is  above  the  steep  slope  or  at  the 
beginning  of  the  gentle  slope.  In  those  two  cases,  the  vertical  directionality  corresponds  to 
a  shallow  source  in  deep  water. 

The  beamforming  results  obtained  with  the  full  water  column  aperture  and  given 
in  Fig.  15,  are  much  more  difficult  to  interpret,  since  actual  arrivals  are  curved  rather  than 
planar  over  the  large  spatial  extent  of  the  vertical  line  array  [14].  The  refraction  index 
causes  an  arrival  to  have  an  angle  varying  locally  as  a  function  of  depth,  according  to 
Snell’s  law.  In  such  conditions,  low  angle  arrivals  can  be  observed,  even  for  a  shallow 
source  in  deep  water,  as  in  the  vertical  directionality  plot  of  the  the  source  above  the 
steep  slope  sector.  The  arrival  structure  is  highly  dependent  on  the  range  between  the 
source  and  the  receiving  array.  When  the  source  is  at  the  beginning  of  the  gentle  slope, 
there  is  no  strong  near  horizontai  arrival.  But,  when  the  source  is  in  shallow  water  or 
above  the  small  2.8  deg  slope,  there  is  significant  power  near  horizontal  in  the  vertical 
directionality  structure,  as  expected. 
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4.  Matched-Field  Processing  Results 

Matched-field  processing  now  is  used  to  study  the  effect  of  downslope  conversion 
and  complement  the  results  presented  in  the  earlier  sections.  The  PE  model  is  used  to 
create  the  set  of  replica  vectors  for  a  search  window  extending  from  25  to  525  km  in  range 
and  from  5  to  985  m  in  depth.  The  MFP  range-depth  cells  are  2  km  wide  and  20  m  deep. 
The  range  step  has  been  chosen  relatively  large  in  order  to  limit  the  number  of  PE  runs.  It 
does  not  need  to  be  significantly  smaller  since  the  interest  in  this  study  is  the  sidelobe 
leakage  into  a  deep  water  MFP  cell  from  shallow  sources  in  the  continental  slope  area. 

For  each  source  range,  PE  is  run  with  an  updated  input  file  to  correctly  account  for 
the  range  dependence  of  the  environment.  The  array  is  assumed  fixed  and  the  motion  of 
the  source  changes  the  source  to  receiver  range  as  well  as  the  environmental  sector  start¬ 
ing  ranges.  Over  the  525  km  range  extent  of  interest,  there  are  ten  sectors  with  different 
sound  speed  profiles.  In  order  to  be  consistent  with  the  pressure  fields  calculated  earlier 
to  simulate  a  source  on  the  continental  slope,  PE  is  started  with  an  initial  field  produced 
by  the  normal  mode  model  KRAKEN,  even  in  the  case  of  deep  sources. 

The  MFP  ambiguity  surfaces  displayed  in  gray-level  correspond  to  the  correlation 
measure  given  by 

Or  in  dB  : 

Pi^s.Zs)  =  10  log  e;.|2  (2) 

where  corresponds  to  the  simulated  data  across  the  vertical  line  array  (due  to  the  shal¬ 
low  source  in  the  continental  slope  area),  normalized  to  unit  norm,  and  C;.  is  the  unit  norm 
replica  vector  which  is  a  function  of  the  hypothesized  source  range  and  depth  and  z^.  H 
denotes  the  Hermitian  (conjugate  transpose)  operator.  This  correlation  measure  actually  is 
a  measure  of  mismatch  between  the  simulated  data  and  the  replica  vectors. 
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4.1.  Processing  the  1000  m  Aperture 

The  processing  results  on  the  1000  m  array  with  80  sensors  are  examined  first. 
The  ambiguity  surface  for  the  source  in  shallow  water,  given  in  Fig.  16,  has  peaks  forming 
patterns  reminiscent  of  a  source  at  depth.  These  features  happen  especially  in  the  275- 
525  km  region.  They  appear  centered  around  a  depth  of  300  m  and  have  a  convergence 
zone  periodicity  on  the  order  of  50  to  60  km.  Hence,  Fig.  16  shows  that  the  shallow  source, 
located  on  the  shelf  in  shallow  water,  1750  km  from  the  receiving  array,  looks  like  a  deep 
source  in  the  MFP  range-depth  search  window.  When  the  source  is  above  the  small  2.8  deg. 
slope,  the  slope  conversion  is  even  stronger.  In  this  case,  the  ambiguity  surface  gray-level 
displayed  in  Fig.  17,  has  the  same  character  as  Fig.  16,  with  even  broader  peaks  distributed 
in  depth  primarily  around  300  m  and  600  m.  A  convergence  zone  periodicity  also  is 
observed.  The  shallow  source  above  the  small  2.8  deg.  slope  does  also  look  like  a  deep 
source  in  the  MFP  search  window. 

The  ambiguity  surfaces  for  the  source  above  the  steep  slopeand  at  the  beginning 
of  the  gentle  slope  are  given  in  Figs.  18  and  19.  Both  surfaces  are  characterized  by  sets  of 
peaks  at  shallow  depths,  around  100  m.  These  surfaces  have  a  character  very  much  dif¬ 
ferent  from  the  two  slope  interacting  source  stations.  They  have  strong  peaks  at  shallower 
depth,  but  none  deeper.  The  shallow  water  sources  in  the  continental  rise  area  (where 
there  is  no  downslope  conversion),  1700  km  from  the  array,  look  like  shallow  sources  in 
the  MFP  search  window. 

To  help  interpret  the  results,  the  beamforming  and  matched-field  processing 
results  are  examined  for  two  sources  located  within  the  MFP  search  window  at  a  range  of 
423  km  and  at  two  different  depths  :  45  m  and  305  m  depth.  The  angular  spectra  obtained 
with  the  1000  m  array  for  the  two  source  depths  are  given  in  Fig.  20.  The  deeper  source 
produces  strong  horizontal  arrivals,  while  the  shallow  source  does  not.  The  MFP  ambiguity 
surfaces  are  gray-level  displayed  in  Figs.  21  and  22.  The  ambiguity  surface  of  the  shallow 
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source  is  characterized  by  small  peaks  distributed  at  shallower  depths  with  a  convergence 
zone  periodicity  in  range.  It  is  similar  to  the  surfaces  where  the  source  is  at  the  beginning 
of  the  gentle  slope  (Fig.  18),  or  in  the  middle  of  the  steep  sector  area  (Fig.  19).  On  the  other 
hand,  the  ambiguity  surface  of  the  deeper  source  produces  butterfly  patterns  with  peaks 
distributed  at  great  depth,  similar  to  the  ambiguity  surfaces  of  the  slope  interacting  sta¬ 
tions. 

The  appearance  of  the  MFP  ambiguity  surface  is  very  different  whether  or  not 
there  is  downslope  conversion,  as  the  initial  transmission  loss  and  conventional  beam- 
forming  results  suggested.  The  shoaling  of  the  sound  axis  seems  to  have  little  effect  on 
the  coupling  of  sound  in  the  deep  sound  channel  for  these  low  frequency  simulations  and 
this  particular  array  configuration.  This  is  consistent  with  the  conventional  beamforming 
results. 

The  location  of  the  highest  peak  in  the  ambiguity  surfaces  and  the  corresponding 
squared  correlation  values  (Eq.  (1)),  are  given  in  Table  1.  The  highest  peak  is  deeper  when 
there  is  downslope  enhancement  than  when  there  is  not.  The  observed  correlation  values 
are  high,  above  0.6  with  the  maximum  value  observed  when  slope  conversion  is  largest  (i.e. 
when  the  source  is  above  the  small  2.8  deg.  slope  sector). 

4.2.  Processing  the  4000  m  Aperture 

A  full  water  column  aperture  array  is  now  assumed  deployed  instead  of  the  1000  m 
array.  This  array  has  323  sensors  with  12.353  m  interelement  spacing.  Initial  beamforming 
results  with  the  large  aperture  array  do  not  yield  dramatic  differences  whether  or  not  the 
source  is  on  the  shelf  slope  or  beyond,  that  is  whether  or  not  there  is  downslope  conver¬ 
sion.  To  resolve  the  structure  of  the  MFP  ambiguity  surface  better,  several  plots  of  ambi¬ 
guity  surfaces  with  different  normalizations  and  dynamic  ranges  are  presented. 

The  ambiguity  surfaces  of  the  four  different  stations,  given  by  Eq.  (2),  are  gray- 
level  displayed  in  Figs.  23  to  26  with  a  20  dB  dynamic  range.  The  two  stations  where  strong 
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slope  conversion  takes  place  (i.e.  when  the  source  is  in  shallow  water  or  above  the  small 
slope),  have  ambiguity  surfaces  with  fairly  broad  peaks.  These  peaks  are  distributed  in 
depth  mostly  around  200  m,  500-600  m  and  800-900m.  As  usual,  convergence  zone  periodi¬ 
city  in  range  is  observed.  When  the  sources  above  the  steep  slope  and  at  the  beginning  of 
the  gentle  slope,  there  is  minimal  downslope  conversion,  and  the  corresponding  MFP  ambi¬ 
guity  surfaces  have  much  greater  granularity.  They  appear  to  have  a  higher  spatial  fre¬ 
quency  content  with  finer  structure  in  both  range  and  depth.  Many  small  peaks  are  distri¬ 
buted  in  range  with  a  convergence  zone  periodicity  which  is  more  pronounced  than  in  the 
case  of  the  slope  interacting  stations.  Also,  the  highest  correlation  levels  appear  to  be 
towards  the  surface  rather  than  larger  depths.  These  observations  are  confirmed  by  the 
plots  of  the  ambiguity  surfaces  normalized  by  their  highest  peak  and  gray-level  displayed 
on  a  10  dB  scale,  in  Figs.  27  to  30. 

As  in  Section  4.1,  it  is  of  interest  to  look  at  the  ambiguity  surfaces  due  to  a  shallow 
(45  m)  and  a  deep  (305  m)  source  within  the  MFP  search  window,  423  km  away  from  the 
receiving  array.  The  ambiguity  surface  due  to  the  shallow  source  is  gray-level  displayed  in 
Fig.  31.  It  is  formed  of  inverted  V  shaped  patterns  with  a  convergence  periodicity  in  range. 
Closer  to  the  array,  the  well  delimited  patterns  are  destroyed  and  tend  to  form  blocks  of 
fine  peaks  peaks  distributed  in  depth  and  range.  As  in  the  case  of  the  1000  m  aperture,  the 
ambiguity  surface  due  to  the  shallow  source  in  the  MFP  search  window  and  the  ambiguity 
surfaces  of  the  non  slope  interacting  stations  have  an  overall  fine  structure  in  common. 
The  ambiguity  surface  when  the  source  is  at  305  m  depth  is  gray-level  displayed  in  Fig.  32. 
The  ambiguity  surface  has  the  familiar  butterfly  pattern  centered  at  the  source  location. 
These  patterns  change  shape  with  range  and  tend  to  spread  at  closer  ranges  to  the  array. 
This  character  is  similar  to  the  one  observed  for  the  slope  interacting  stations  (see  Figs.  29 
and  30). 


Reducing  the  dynamic  range  of  the  gray-level  plot  from  20  down  to  10  dB  has  a 
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clipping  effect  and  removes  most  of  the  lower  correlation  level  structure  in  the  surface 
given  by  Eq.(l).  These  ambiguity  surfaces  given  in  Figs.  33  to  36  give  an  opportunity  of 
detecting  the  main  structural  changes  in  the  ambiguity  fields  when  the  source  does  or 
does  not  interact  with  the  slope.  Indeed,  one  can  observe  in  the  case  of  the  source  in  shal¬ 
low  water  and  above  the  small  slope,  peaks  distributed  in  depth  around  200  m  as  well  as 
deeper,  below  500  m.  These  deep  ambiguity  peaks  are  not  observed  in  the  ambiguity  sur¬ 
faces  of  the  sources  above  the  steep  slope  and  at  the  beginning  of  the  gentle  slope  sector. 
In  those  surfaces,  only  a  few  peaks  around  100  m  are  observed.  It  is  of  interest  to  com¬ 
pare  these  surfaces  with  the  ones  (with  the  same  normalization)  corresponding  to  a  source 
within  the  range-depth  MFP  window,  at  423  km  range  at  two  different  depths  :  45  m  or  305 
m  depth.  The  ambiguity  surfaces  are  given  in  Figs.  37  and  38.  There  is  a  strong  similarity 
between  the  ambiguity  surface  of  the  source  at  45  m  depth  and  423  km  range  and  the  one 
of  the  non-slope  interacting  stations.  There  is  as  well  a  strong  similarity  between  the  ambi¬ 
guity  surface  of  the  source  at  305  m  depth  and  423  km  range  and  the  one  of  the  slope 
Interacting  stations. 

Table  2  provides  quantitative  results  of  matched-field  processing  on  the  full  water 
column  array.  The  correlation  values  for  this  323  sensor  array  are  still  high,  around  0.3. 
The  depth  of  the  highest  peak  does  not  correlate  as  clearly  with  the  slope  interaction,  as  in 
the  case  of  the  1000  m  aperture.  This  is  consistent  with  the  fact  that  the  ambiguity  sur¬ 
faces  are  not  dramatically  different  whether  or  not  there  is  downslope  conversion  as  they 
are  for  the  1000  m  aperture. 

5.  Conclusions 

The  manifestation  of  sources  on  the  continental  slope  in  the  deep  water  MFP  ambi¬ 
guity  surface  has  been  studied  in  this  report.  The  PE  model  is  used  to  simulate  the  pres¬ 
sure  fields  due  to  a  low  frequency  (18  Hz)  shallow  source  that  has  propagated  over  1750 
km.  Two  receiving  arrays  are  considered;  (1)  a  1000  m  aperture  with  80  sensors  deployed 
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from  400  to  1390  m  and  (2)  a  full  water  column  array  with  323  sensors.  In  both  cases,  the 
interelement  spacing  is  constant  and  equal  to  12.353  m.  The  environment  is  range  depen¬ 
dent  in  both  bathymetry  and  sound  speed.  After  a  preliminary  analysis  of  the  simulated 
pressure  fields  was  carried  out  using  transmission  loss  plots,  field  plots  and  conventional 
plane-wave  beamforming,  matched-field  ambiguity  surfaces  were  calculated  on  a  range- 
depth  search  window  extending  in  range  from  25  to  525  km  and  in  depth  from  5  to  985  m. 
These  surfaces  were  studied  when  the  source  was  (1)  on  the  continental  shelf  in  shallow 
water,  1750  km  from  the  receiving  array,  (2)  above  a  small  2.8  deg.  slope,  (3)  above  a  steep 
1 1.5  deg.  slope,  and  finally  (4)  in  deep  water  at  the  beginning  of  a  gentle  0.4  deg.  slope  sec¬ 
tor  where  the  deep  water  begins. 

Array  design  has  an  impact  on  the  appearance  of  the  MFP  ambiguity  surfaces.  The 
shorter  array  (1000  m)  deployed  across  the  sound  axis  yields  very  different  results 
whether  or  not  the  source  is  at  a  location  where  downslope  conversion  takes  place.  When 
the  source  is  above  the  small  slope  sector  and  in  shallow  water,  the  MFP  ambiguity  sur¬ 
faces  have  broad  peaks  at  large  depths.  The  shallow  source  looks  like  a  deep  source 
because  of  downslope  conversion.  On  the  other  hand,  when  the  source  is  above  the  steep 
slope  or  at  the  beginning  of  the  gentle  slope,  there  is  no  slope  interaction  and  the  ambi¬ 
guity  surfaces  have  smaller  peaks  at  shallow  depth.  The  shallow  source  still  looks  like  a 
shallow  source.  For  this  array  configuration,  both  conventional  beamforming  and 
matched-field  processing  results  show  that  shoaling  of  the  sound  axis  towards  the  con¬ 
tinental  slope  (which  is  another  proposed  mechanism  for  the  coupling  of  sound  in  the 
deep  sound  channel  [5,8])  has  little  impact. 

The  results  obtained  with  the  4000  m  array  (full  water  column  aperture)  are  not  as 
dramatic  as  the  shorter  array  results,  although  a  careful  analysis  of  the  MFP  ambiguity  sur¬ 
faces  lead  to  the  observation  of  definite  differences  between  the  slope  interacting  and 
non-interacting  sources.  The  character  of  the  MFP  ambiguity  surfaces  for  sources  which  do 
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not  interact  with  the  slope  and  of  a  shallow  source  in  deep  water  are  similar.  They  have  a 
fine  structure  and  a  concentration  of  the  strongest  peaks  towards  the  surface.  On  the 
other  hand,  the  ambiguity  surfaces  of  the  slope  interacting  sources  have  larger  features 
distributed  at  larger  depths.  These  latter  are  similar  to  the  ambiguity  surfaces  of  deep 
sources  in  deep  water. 

For  both  array  apertures,  high  MFP  correlation  values  (Eq.  (1))  are  observed  with 
values  between  0.6  and  0.7  for  the  1000  m  array  and  between  0.23  and  0.32  for  the  full 
water  column  array.  This  means  that  an  18  Hz  source  on  the  continental  shelf  and  at  large 
range  can  look  like  a  source  in  deep  water  at  closer  range  to  the  array.  The  highest  corre¬ 
lation  values  are  observed  when  the  shallow  source  is  above  the  small  2.8  deg.  slope  where 
there  is  maximum  slope  interaction. 
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Table  1.  Quantitative  results  of  matched-field  processing  on  the  1000  m  array. 


Station 

Correlation 

Depth  (m) 

Range  (km) 

Shallow  water 

0.64303 

565 

513 

Small  Slope 

0.69464 

725 

525 

Steep  Slope 

0.66555 

105 

351 

Gentle  Slope 

0.65828 

105 

203 

Table  2.  Quantitative  results  of  matched-field  processing  on  the  4000  m  array. 


Station 

Correlation 

Depth  (m) 

Range  (km) 

Shallow  water 

0.64303 

565 

513 

Small  Slope 

0.69464 

725 

525 

Steep  Slope 

0.66555 

105 

351 

Gentle  Slope 

0.65828 

105 

203 
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List  of  Figures 

Fig.  1  :  Overall  propagation  geometry  of  the  simulation. 

Fig.  2  :  Sound  speed  field  in  deep  water.  The  receiving  array  is  on  the  left  and  the  con¬ 
tinental  slope  on  the  right.  The  sound  speed  profiles  of  the  28  environmental  sectors  are 
shifted  by  5  m/s  per  sector. 

Fig.  3  ;  Sound  speed  field  and  bathymetry  in  the  slope  area. 

Fig.  4  :  Bottom  sound  speed  profiles  in  the  slope  area. 

Fig.  5  :  Transmission  loss  versus  range  for  the  four  stations  on  the  slope.  The  source  is  45 
m  deep  and  the  receiver  is  1000  m  deep.  The  receiving  array  is  at  range  1750  km.  The  top 
curve  corresponds  to  the  source  above  the  small  2.8  deg  slope. 

Fig.  6  :  Field  plot  for  the  45  m  deep  source  on  the  shelf. 

Fig.  7  :  Bathymetry  corresponding  to  Fig.  6.  The  source  is  on  the  shelf. 

Fig.  8  :  Field  plot  for  the  45  m  deep  source  above  the  small  2.8  deg.  slope. 

Fig.  9  ;  Bathymetry  corresponding  to  Fig.  8.  The  source  is  above  the  small  slope. 

Fig.  10  :  Field  plot  for  the  45  m  deep  source  above  the  steep  1 1.5  deg.  slope  sector. 

Fig.  11  :  Bathymetry  corresponding  to  Fig.  10.  The  source  is  above  the  steep  slope. 

Fig.  12  :  Field  plot  for  the  45  m  deep  source  in  deep  water. 

Fig.  13  ;  Bathymetry  corresponding  to  Fig.  12.  The  source  is  in  deep  water. 

Fig.  14  :  Angular  spectra  obtained  with  conventional  beamforming  on  the  1000  m  aperture. 
The  45  m  deep  source  is  located  at  the  four  different  locations  on  the  slope. 

Fig.  15  :  Angular  spectra  obtained  with  conventional  beamforming  on  the  4000  m  aperture. 
The  45  m  deep  source  is  located  at  the  four  different  locations  on  the  slope. 

Fig.  16  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  on  the  shelf,  1750  km 
from  the  1000  m  long  array. 

Fig.  17  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  on  the  small  2.8  deg. 
slope,  1729.6  km  from  the  1000  m  long  array. 

Fig.  18  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  on  the  the  steep  11.5  deg. 
slope,  171 1.9  km  from  the  1000  m  long  array. 

Fig.  19  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  at  the  the  beginning  of 
the  gentle  0.4  deg.  slope,  1707.5  km  from  the  1000  m  long  array. 

Fig.  20  :  Angular  spectra  obtained  with  conventional  beamforming  on  the  1000  m  aperture. 
The  source  is  located  within  the  MFP  range-depth  search  window  at  423  km  range  and  at 
two  different  depths  :  45  m  and  305  m. 

Fig.  21  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  423  km  from  the  1000  m 
long  array. 

Fig.  22  :  MFP  ambiguity  surface  for  the  305  m  deep  source  located  423  km  from  the  1000  m 
long  array. 

Fig.  23  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  on  the  continental  shelf, 
1750  km  from  the  4000  m  long  array. 

Fig.  24  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  on  the  the  small  2.8  deg. 
slope,  1729.6  km  from  the  4000  m  long  array. 

Fig.  25  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  on  the  steep  11.5  deg. 
slope,  1711.9  km  from  the  4000  m  long  array. 

Fig.  26  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  at  the  beginning  of  the 
gentle  0.4  deg.  slope,  1707.5  km  from  the  4000  m  long  array. 
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Fig.  27  :  Normalized  MFP  ambiguity  surface  for  the  45  m  deep  source  located  on  the  con¬ 
tinental  shelf,  1750  km  from  the  4000  m  long  array. 

Fig.  28  :  Normalized  MFP  ambiguity  surface  for  the  45  m  deep  source  located  on  the  small 
2.8  deg.  slope,  1729.6  km  from  the  4000  m  long  array. 

Fig.  29  :  Normalized  MFP  ambiguity  surface  for  the  45  m  deep  source  located  on  the  steep 
1 1.5  deg.  slope,  171 1.9  km  from  the  4000  m  long  array. 

Fig.  30  :  Normalized  MFP  ambiguity  surface  for  the  45  m  deep  source  located  at  the  begin¬ 
ning  of  the  gentle  0.4  deg.  slope,  1707.5  km  from  the  4000  m  long  array. 

Fig.  31  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  423  km  from  the  4000  m 
long  array. 

Fig.  32  :  MFP  ambiguity  surface  for  the  305  m  deep  source  located  423  km  from  the  4000  m 
long  array. 

Fig.  33  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  on  the  continental  shelf, 
1750  km  from  the  4000  m  long  array.  The  dynamic  range  is  reduced  to  10  dB. 

Fig.  34  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  on  the  small  2.8  deg. 
slope,  1729.6  km  from  the  4000  m  array.  The  dynamic  range  is  reduced  to  10  dB. 

Fig.  35  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  on  the  the  steep  11.5  deg. 
slope,  1711.9  km  from  the  4000  m  long  array.  The  dynamic  range  is  reduced  to  10  dB. 

Fig.  36  :  MFP  ambiguity  surface  for  the  45  m  deep  source  located  at  the  beginning  of  the 
gentle  0.4  deg.  slope,  1707.5  km  from  the  4000  m  long  array.  The  dynamic  range  is  reduced 
to  10  dB. 

Fig.  37  ;  MFP  ambiguity  surface  for  the  45  m  deep  source  located  423  km  from  the  4000  m 
long  array.  The  dynamic  range  is  reduced  to  10  dB. 

Fig.  38  :  MFP  ambiguity  surface  for  the  305  m  deep  source  located  423  km  from  the  4000  m 
long  array.  The  dynamic  range  is  reduced  to  10  dB. 
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Bathymeiry  :  Source  in  small  slope  secior 


0  25  50  75  100  125  150 


Range  (Un) 


Depth  (10  *  m) 


Figure  11 


Bathymetry  :  Source  in  steep  slope  sector 
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Bathymeiry  :  Source  in  genile  slope  sector 
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Figure  17 

Source  out  of  the  MFP  window  at  45  m  depth 
1000  m  array  with  12.353  m  interelement  spacing 
Small  slope  sector 
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Figure  18 

Source  out  of  the  MFP  window  at  45  m  depth 
1000  m  array  with  12.353  m  interelement  spacing 
Steep  slope  sector 
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Figure  19 

Source  out  of  the  MFP  window  at  45  m  depth 
1000  m  array  with  12.353  m  interelement  spacing 
Gentle  slope  sector 
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Source  oi  45  m  depth  and  423  kn  range 

1000  m  array  with  12.353  m  interelement  spacing 
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Figure  22 


Source  at  305  m  depth  and  423  kn  range 

1000  m  array  with  12.353  m  interelement  spacing 


Source  out  of  the  MFP  window  at  45  m  depth 
4000  m  array  with  12.353  m  interelement  spacing 
Shallow  water  sector 
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Source  out  of  the  MFP  window  at  45  m  depth 
4000  m  array  with  12.353  m  interelement  spacing 
Small  slope  sector 
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Figure  25 

Source  out  of  the  MFP  window  at  45  m  depth 
4000  m  array  with  12.353  m  interelement  spacing 
Steep  slope  sector 
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Source  out  of  the  MFP  window  at  45  m  depth 
4000  m  array  with  12.353  m  inter  element  spacing 
Gentle  slope  sector 
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Figure  27 


Source  out  of  the  MFP  window  at  45  m  depth 
4000  m  array  with  12.353  m  interelement  spacing 
Shallow  water  sector 
Surface  normalized  to  highest  peak 


Source  out  of  the  MFP  window  at  45  m  depth 
4000  m  array  with  12.353  m  interelement  spacing 
Small  slope  sector 
Surface  normalized  to  highest  peak 
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Figure  29 


Source  out  of  the  MFP  window  at  45  m  depth 
4000  m  array  with  12.353  m  interelement  spacing 
Steep  slape  sector 
Surface  normalized  to  highest  peak 
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Source  out  of  the  MFP  window  at  45  m  depth 
4000  m  orray  with  12.353  m  interelement  spacing 
Gentle  slope  sector 
Surface  normalized  to  highest  peak 


275  300  325  350  375  400  425  450  475  500  525 

Range  (Wn) 


-10  -6  -6  -4  -2  0 


Relative  Power  (dB) 


Depth  (10  *  m)  Depth  (10 


Figure  31 


Source  ai  45  m  depth  and  423  Un  range 

4000  m  array  with  12.353  m  interelement  spacing 
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Source  oi  305  m  depth  and  423  Un  range 

4000  m  array  with  12.353  m  interelement  spacing 
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Source  out  of  the  MFP  window  at  45  m  depth 
4000  m  array  with  12.353  m  interelement  spacing 
Shallow  water  sector 
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Source  out  of  the  MFP  window  at  45  m  depth 
4000  m  array  with  12.353  m  interelement  spacing 
Small  slope  sector 
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Figure  35 

Source  out  of  the  MFP  window  at  45  m  depth 
4000  m  array  with  12.353  m  interelement  spacing 
Steep  slope  sector 
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Figure  36 

Source  out  of  the  MFP  window  at  45  m  depth 
4000  m  array  with  12.353  m  interelement  spacing 
Gentle  slope  sector 
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Figure  37 

Source  at  45  m  depth  and  423  kn  range 

4000  m  array  with  12.353  m  interelement  spacing 


Figure  38 

Source  oi  305  m  depth  and  423  kn  range 

4000  m  array  with  12.353  m  interelement  spacing 
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